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We have developed a method to equip homodyne interferometers with the capability to operate 
with constant high sensitivity over many fringes for continuous real-time tracking. The method 
can be considered as an extension of the " Ji . . . J4" methods, and its enhancement to deliver very 
sensitive angular measurements through Differential Wavefront Sensing is straightforward. Beam 
generation requires a sinusoidal phase modulation of several radians in one interferometer arm. On 
a stable optical bench, we have demonstrated a long-term sensitivity over thousands of seconds of 
ClmradZ-s/Hz that correspond to 20pm/v'Hz in length, and 10nrad/\/Hz in angle at millihertz 
frequencies. 
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I. INTRODUCTION 

Optical interferometers with sub-wavelength resolu- 
tion are useful in many optical metrology applications, 
such as, for example, length measurements, gravitational 
wave detection, wavefront sensing, and surface profiling, 
among others. Our technique was developed in the con- 
text of continuously measuring the position and orienta- 
tion of a free-floating test mass for space-based gravita- 
tional wave detection [T] , although the method is useful 
for other applications as well. Other techniques for the 
optical readout of free-floating test masses at millihertz 
frequencies are currently under investigation, such as a 
polarizing heterodyne interferometer reaching a sensitiv- 
ity of about 300pm/vHz [5], a compact homodyne in- 
terferometer with a sensitivity of 100pm/'\/Hz [3], and a 
robust implementation of an optical lever with a read- 
out noise level of 100pm/vHz[T. Another method to 
do this is heterodyne interferometry as developed for 
LISA Pathflnder [S] with a sensitivity of better than 
5pm/vHz'S'. The method we present here achieves an 
optical pathlength measurement sensitivity of the order 
of 20 pm/vHz, and with an angular resolution better 
than lOnrad/vHz, both above 3mHz. The conversion 
from real test mass motion to optical pathlength is given 
by the interferometer topology, and is in our case about 
a factor of 2, which yields a test mass motion resolu- 
tion of approximately lOpm/vHz. Those interferome- 
ters with the highest accuracy, namely Fabry-Perot in- 
terferometers on resonance or recycled Michelson inter- 
ferometers on a dark fringe [7j , have a dynamic range 
of a small fraction of one fringe only. High resolution 
and wide dynamic range can be simultaneously achieved 
by, e.g., active feedback or heterodyning, each of which 
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has disadvantages. Active feedback transfers the inher- 
ent non-linearity of the feedback actuator to the output 
signal or requires another stabilized laser and a measure- 
ment of the high-frequency beat note. Heterodyning, 
on the other hand, requires a complex setup to gener- 
ate the two coherent beams with a constant frequency 
difference, typically involving two acousto-optic modula- 
tors (AOMs) with associated frequency generation and 
RF power amplification. Other methods to overcome 
these limitations involve variations of sinusoidal phase 
shifting interferometry [8l-[T2] , reporting accuracies of the 
order of 1 nm. These methods are typically used in 
"single-shot" mode for static applications such as sur- 
face profiling, whereas our method is designed for con- 
tinuous real-time, long-term tracking of a moving target 
with low noise at millihertz frequencies. In particular, 
the so-called "Ji . . . J4" method [TBHT^ . involves a sinu- 
soidal phase modulation at a fixed frequency /mod with 
modulation depths m « 1 ... 5 in one arm of the inter- 
ferometer. The spectrum of the resulting photocurrent 
has components at integer multiples of /mod, with am- 
plitudes that can be written in terms of the Bessel func- 
tions J„(to) (hence the name) and the phase difference 
ip due to the optical pathlength difference. The meth- 
ods then proceed to use analytical formulae to solve for 
the unknowns m and ip, after obtaining the harmonic 
amplitudes from a spectrum analyzer or a Fast Fourier 
Transform (FFT) of the digitized time series. The accu- 
racies reported are of order 10. . . 100 mrad (1.7 .. . 17 nm) 
for a laser wavelength of 1064 nm. We generalize this ap- 
proach by using a higher modulation index m (up to 10 
or 20) and making use of all harmonics up to an order 
N fn m. These are more observations than the four un- 
knowns (m, ip, modulation phase and a common factor), 
making an analytical solution impossible. Instead we use 
a numerical least-squares solution which allows consis- 
tency checks and improves the signal-to-noise ratio. For 
our typical applications we keep m near constant at an 
optimal value and take ip as useful output, achieving an 
accuracy of better than 0.12mrad/\/Hz that correspond 



to 20pm/vHz in length, and lOnrad/vHz in angle at 
millihertz frequencies. As compared to heterodyne inter- 
ferometers, more complex data processing is necessary to 
recover the optical pathlength from the measured pho- 
tocurrent. However, with the availability of inexpensive 
processing power, this computational complexity is often 
preferable to additional optics and electronics hardware 
needed for the optical heterodyning. 



II. THEORY 

The signal Vp-oit) of a photodetector at the output 
of a phase-modulated homodyne interferometer can be 
expressed as 

Vpoit) = A [1 + C cos {tp + m cos {uj,-nt + ip))] , (1) 

where ip is the interferometer phase, m is the modulation 
depth, Win = 27r /m is the modulation frequency, ^p is 
the modulation phase, C < 1 is the contrast, and A com- 
bines nominally constant factors such as light powers and 
photodiode efficiencies. The interferometer output is pe- 
riodic with /m and its signal waveform characteristically 
depends on the interferometer phase (p. Figure [l] illus- 
trates typical waveforms obtained for various states of p. 
The expression of Equation [l] can be expanded into its 
harmonic components as: 



VpD{t) = Vdc(<p) + ^a„(TO, (p)cos(n(a;mt + V')) (2) 



n=l 



with 



an{m,(p) — k Jn{m) cos iip + n—j , and 

Vdc('/') = a {1 + CJoim) cos If), 



(3) 
(4) 



where k = 2CA, and Jnifn) arc the Bcsscl functions. 
Figure [2] shows the dependence of the harmonic ampli- 
tudes an{m, if) in terms of (p. Our technique is centered 
around these harmonic amplitudes anim, ip) which on the 
one hand can be directly measured by numerical Fourier 
analysis of the photocurrent, and on the other hand have 
the above analytical relationships to the unknowns p, m, 
-tp, k. The technique we present here uses higher modula- 
tion depths 771 > 6 to set up an overdimensioned system 
of equations that can be numerically solved for the four 
sought parameters p, m, ip, and the common factor k 
by a least-squares fit algorithm. The information of the 
harmonic amplitude aQ{m,ip), corresponding to the DC 
component Vaci'P) i^ ^ot used by the fit algorithm, since 
it usually contains a higher noise level due to large vari- 
ations in environmental and equipment conditions such 
as room illumination and electronic noise, among others. 
However, it is useful for computation of the interferome- 
ter visibility and alignment signals. 



III. DATA PROCESSING 

The signal Vpd(0 measured at the photodetector is 
digitized after appropriate analog processing and anti- 
alias filtering. The sampling rate /samp is arranged to be 
coherent to the modulation frequency /mod- The time 
series is split in segments of length iVppT samples that 
are processed by a Discrete Fourier Transform in order to 
compute N = N-pyt/2 complex amplitudes q;„(to, 93). A 
non-linear fit algorithm is applied to match the measured 
an{m,p)) to the complex amplitudes c„ computed from 
the model 



a„(m, p) = anim, p) e™'^ . 



(5) 



There is a total of 2N equations that can be set up in 
two uncorrelated system of equations: 



nip — arctan 



/3{a„(m, v?)} 



VK{a„(m, 93)} 
a„{m, p) = a„(m, p) e"'"'^, n = 1, 2, 3 



n = 1, 2, 3 . . . iV, (6) 

N, (7) 



where an{m,p) e^™"^ is a real number. For the mea- 
sured a„(m, (/j) e"'"''', this is not exactly the case due 
to noise and phase distortions introduced by the analog 
electronics. In order to solve the system of equations, a 
Levenberg-Marquardt fit algorithm (THl [IT] is applied to 
minimize the least-squares expression 



N 



y^ (an("i, p) - an{m, p)Y 



(8) 



where x^ is a four dimensional function of to, p, ip, and k. 
In practice, these parameters barely vary between con- 
secutive segments of length TVfft, giving good starting 
values for a rapid convergence of the fit. Only in the case 
this is not accomplished such as upon initialization or 
after large disturbances, a modified version of the more 
robust Nelder-Mead Simplex algorithm [T5] is applied as 
initial step. In order to find best values of the modu- 
lation index m and the number of bins N for optimum 
performance, we conducted a numerical analysis of the 
4x4 Hessian matrix of x^ that is given by 



H = (i?,,) 



2^,2 



d^X 



dfljd^o 



(9) 



where fi = {to, p, ip, k} are the four parameters. The 
inverse of the Hessian matrix H~^ — (rjij) yields infor- 
mation about the parameter estimates on the variances 



a and correlation coefficients 



Pij- 






Pi3 



(10) 
(11) 



An excursion of p over the range [0, 27r] -which corre- 
sponds to one interferometer fringe- was conducted in 
64 steps by fixed N and m in order to compute the 
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FIG. 1. Waveform of the obtained interferogram for different operating points of the interferometer phase ^ with a modulation 
depth m = 6 rad. 
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FIG. 2. Dependence of the harmonics ampUtudes a„{m,ip) 
with respect to the interferometer phase if with a modulation 
depth 771 = 6 rad. 
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FIG. 3. Ideal resolution in (fi as function of the modulation 
index m for A'^ = 10, for the best and worst ip as well as the 
average for all (p G [0, 2n]. 



best, worst and average values of the standard deviation 
an-lNjirijip), which are shown in Figure pj Assuming 
worst case values of the variances 

o'o,(^:H= max o-^^(Ar,TO,(/?), (12) 

ipe[0,2ir] 

we run a similar analysis varying N and m to evaluate 
for best resolution of any value of (/?, which is our main 
measurement and often not entirely under control. The 
results are shown in Figure |4] This analysis revealed 
useful parameter estimates for 3 < m < N, and pos- 



sible best values of m for minimum (T^p in the cases of 
?7i = 6, A^ > 8 and m = 9,N > 10, suggesting best res- 
olutions of (p. These results are only rough guidelines, 
since real instrument noise has not been yet considered. 
The dominant noise sources have, however, been investi- 
gated experimentally, as discussed in S ect ion |V] below. In 
addition, software simulations of the fit routine were run 
with synthetic data as input. Hardware characteristics of 
the data acquisition system (DAQ) such as digitization 
effects and frequency response of the anti-aliasing filter 
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FIG. 4. Ideal resolution in ip as function of the modulation 
index m for different orders A'^, for the worst value of tp at 
each point of each curve. 



were considered in the generation of mock-data, using 
Equation [l] as nominal noise- free model. We introduced 
two independent mock-data sets into two virtual DAQ 
channels, by linearly increasing m with A^ = 10 bins, 
and recorded the fit output (phase Lp). We then com- 
puted their phase difference and extracted the nominal 
offset in oder to obtain the dependence of the phase fluc- 
tuations (noise) with ttt,, which is shown in Figure [5] A 
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FIG. 5. Dependence of the measured phase noise with m. 



minimum can be observed around m = 9.5 — 10, which is 
consistent with the analysis presented in Figure]!] Hence, 
a DAQ test system for real optical length measurements 
was set up with A'^ = 10 bins, and a modulation index 
m«9.7. 



IV. EXPERIMENTAL SETUP 

We have applied this technique to a very stable in- 
terferometer, namely the engineering model of the LISA 
Pathfinder (LPF) optical bench [5] , which consists of a 
20cmx20cm Zerodur® baseplate with optical compo- 
nents fixed by hydroxide-catalysis bonding |19] . This op- 
tical bench has been extensively characterized as part of 
ground testing campaigns for the optical metrology of 
the LISA Pathfinder mission [5D] , and its optical path- 
length stability has been measured to be better than 
5pm/vHz above 1 mHz. A non-planar ring oscillator 
(NPRO) Nd:YAG laser producing 300 mW at 1064 nm 
was used as light source. For the experimental test, 
we chose a two-beam Mach-Zehnder interferometer, us- 
ing self-assembled fiber-coupled phase modulators con- 
sisting of single-mode fiber optics coiled around ring 
piezo-electric transducers (RPZT) in order to reach high 
modulation depths (up to 10 or 20). Figure |6] shows a 
schematic overview of the setup. The laser beam is split 
into two equal parts at the first beamsplitter BSl. A 
RPZT driven by a sinusoidal voltage of approximately 
4.5 Vpp at /mod = 280 Hz, produces a phase modula- 
tion of modulation depth m « 9.7 in one of the two 
beams. This portion of the optical setup denoted as 
modulation bench, contains the first beamsplitter BSl, 
phase modulator, and corresponding fiber coupling de- 
vices which are all mounted on a standard metal optical 
breadboard. A single-mode fiber feed-through is used to 
bring the main laser beam into a vacuum chamber where 
both, the modulation bench and the optical bench reside. 
The LISA Pathfinder optical bench is a set of four non- 
polarizing Mach-Zehnder interferometers, three of which 
have been used in these experiments. The first one - 
denoted M- measures distance fluctuations between two 
mirrors mounted on 3-axes piezo-electric actuators [6J . A 
second one -denoted R- serves as phase reference to 
cancel common-mode pathlength fluctuations that arise 
at the modulation bench, such as in metal mounts, phase 
modulator, and fiber optics. The third interferometer - 
denoted F- has an intentionally large optical pathlength 
difference of approximately 38 cm, and is used to measure 
laser frequency fiuctuations. If we denote by sm and 
sr the optical pathlengths of the measurement and ref- 
erence interferometer, respectively, the phases emerging 
from the fit algorithm are given by 

VM = ¥ {(Sl + sm) - (S2 + S3)} - ^ {{sM - S3) + A} (13) 
VR^^ {(Sl + Sr) - {S2 + Si)} = ^ {{SR - S4) + A} ,(14) 

where A = 1064 nm is the laser wavelength, s^ are the 
optical paths outlined in Figure ]6] and A = si — S2 rep- 
resents the common-mode pathlength difference between 
the two beams that includes everything starting from the 
first beamsplitter BSl, the modulator, fiber optics up to 
the beamsplitters on the stable optical bench. Typically, 
the fiuctuations of A are several /im on 10. . . 1000 second 
time scales and thus much larger than what we want to 
measure. However, the optical pathlengths sr, S3 and S4 
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FIG. 6. Schematic overview of the experiment. 



are confined to the stable optical bench and have only 
negligible fluctuations. By measuring both ipM and (pn 
and computing their difference 



<(2 = ^AI ~ ^R 



27r 



{SM - {sr + S3~ S4)} (15) 



it is possible to cancel the common-mode fluctuations A 
and to obtain a measurement that is dominated by the 
fluctuations of sm as desired. All photodetectors are in- 
dium gallium arsenide (InGaAs) quadrant diodes with 
5 mm diameter. The photocurrent of each quadrant is 
converted to a voltage with a low-noise transimpedance 
amplifier, filtered with a 9-pole 8 kHz Tschebyscheff anti- 
aliasing filter and digitized at a rate /samp — 20 kHz by a 
commercial 16-channel, 16-bit analog-to-digital converter 
(ADC) card installed in a standard PC running Linux. 
The time series are split in segments of TVpFT = 1000 
samples and transformed by a FFT algorithm [3T] . The 
iV = 10 complex amplitudes of bins 1. . . 10 of /mod at 
frequencies 280 . . . 2800 Hz are then fitted. This configu- 
ration allows us to reach a real-time phase measurement 
rate /^ = /samp/A^FFT = 20 Hz. 



V. NOISE INVESTIGATIONS 

During test and debugging experiments, two main 
noise sources were identified to limit the interferometer 
sensitivity with this technique, which are laser frequency 
noise, and the frequency response (transfer function) of 
the DAQ analog electronics, including photodiode tran- 
simpedance amplifiers and anti-aliasing filters. In the fol- 
lowing we explain the coupling mechanism of these noise 
sources, and the mitigation strategies we implemented to 
counteract them. 



A. Laser frequency noise 

Laser frequency noise translates into phase readout 
noise in any interferometer, whose pathlength difference 
As between the two interfering beams is not exactly 
zero. In the case of the LPF optical bench, this path- 
length mismatch has been determined to be approxi- 
mately 10mm[5j. The free-running frequency noise Sv 
of an unstabilized Nd:YAG NPRO laser at 10 mHz has 
been measured to be of the order of 2 x 10^ Hz/\/lIz '22|. 
The conversion factor from laser frequency fluctuations 
dv into phase fluctuations Sip is given by the difference 
in time of travel between the two beams A.s/c, such that 
an estimate of the noise level can be calculated as 



^ c 3 X 108 m/s 



2 X 10^ Hz = 0.4 mrad/ VHz, 
(16) 



which limits the interferometer optical pathlength reso- 
lution 5s to 



A 



6s = — 5w = 
27r ^ 



1064 nm 

2^ 



0.4mrad/vHz = 68pm/vHz. 

(17) 

We implemented two mitigations strategies to correct for 
this effect and improve the length resolution. Both meth- 
ods worked similarly well and allowed suppression of this 
error below the other noise terms. The first one is based 
on the active laser frequency stabilization, for which 
we have used a commercial iodine-stabilized Nd:YAG 
laser. The second method uses the third interferometer 
F (mentioned above) to independently measure a phase 
proportional to the amplified laser frequency fluctua- 
tions, and applies a noise subtraction technique |6ll23t[24] 
that properly estimates the coupling factor and removes 
the contribution from the final data stream. The phase 



of this interferometer was read out with the same deep 
modulation method as the main channels, and is dom- 
inated by laser frequency fluctuations due to its large 
pathlength difference (« 38 cm). A third method can 
also be easily implemented as an active stabilization loop 
by feeding back to the laser over a digital-to-analog con- 
verter the output of a digital controller that uses the 
difference phase extracted from interferometers F and R 
as error signal. 



B. Frequency response of data acquisition system 

The dominant error was identified to be the frequency 
response of the analog electronics of the data acquisition 
system, in particular the contribution of the photodiode 
transimpedance amplifiers and the anti-aliasing filters. 
The transfer function (TF) of this analog portion of the 
DAQ shows small ripples in its magnitude of the order 
of 0.9 dB. The desired parameter (p is essentially deter- 
mined by running a fit onto the relative amplitudes of 
the 10 measured harmonic components an{m^Lp). The 
ripples are, however, large enough to alter the ratio be- 
tween the harmonic amplitudes, such that the fit algo- 
rithm is disturbed, resulting in a high noise level. We 
removed this error by separately measuring the transfer 
function of each channel of the analog front end, fitting 
it to a model, and correcting accordingly the measured 
complex amplitudes an{m^ if) before entering the fit rou- 
tine. Thus, we obtained the corresponding TF complex 
values Pn — bn e*^" for the 10 frequency bins of interest 
280 — 2800 Hz. Hence, the measured complex amplitudes 
an{'m,(p) were corrected as 



a'^{m,ip) 



an{m,ip) 



(18) 



By using complex numbers, this correction also accounts 
for the TF phase shift, and improves the estimation ca- 
pability of the modulation phase ip- 



VI. OPTICAL LENGTH AND ATTITUDE 
MEASUREMENTS 
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FIG. 7. Sensitivity of real optical pathlength measurements. 
Dashed curve with crosses: initial sensitivity prior to noise 
correction techniques. Dashed curve: sensitivity upon cor- 
rection of DAQ frequency response. Solid curve: sensitivity 
reach after application of noise mitigation strategies -laser fre- 
quency noise and DAQ frequency response-. 



an additional factor of approximately 3.5 at 10 niHz. The 
measured optical pathlength sensitivity of this technique 
is of the order of 20pm/vHz above 3mHz and approx- 
imately a factor of 2 above the performance required to 
the LPF interferometry, which has been plotted for com- 
parison purposes. As mentioned above, all photodetec- 
tors at the interferometer outputs on the LPF optical 
bench are quadrant cell diodes. The phases extracted 
from each individual quadrant cell are processed by a 
differential wavefront sensing (DWS) algorithm |251 [26]. 
in order to measure the interferometer alignment with 
high angular resolution. The results of this measurement 
are shown on Figure |8] as a linear spectral density. As it 
can be read from the plot, this technique reaches an an- 
gular sensitivity better than lOnrad/vHz above 3mHz, 
meeting with sufficient margin the requirements set to 
the LPF interferometry that have been also included in 
the graph as a comparison. 



The experimental setup of Figure [6] was used to con- 
duct long-term interferometric length measurements on 
the LPF optical bench. Figure [7] shows the results ob- 
tained in form of linear spectral densities. The dashed 
curve with crosses is the sensitivity obtained initially with 
this method, without applying any of the noise mitigation 
strategies explained in Section |V] The dashed curve is 
the sensitivity achieved after applying the complex value 
correction of the DAQ frequency response to the mea- 
sured harmonic amplitudes a„(TO,(/5) (as given by Equa- 
tion 18), resulting in a sensitivity improvement of about 
one order of magnitude. The solid curve is the measure- 
ment length sensitivity reached upon subtraction of laser 
frequency noise, which increases the length resolution in 



VII. COMPARISON WITH OTHER 
TECHNIQUES 

The only method known to the authors that allows 
length and angular measurements at arbitrary operating 
points with low noise at millihertz frequencies is het- 
erodyne interferometry as described in Ref. [S]. The 
deep phase modulation method presented here needs, 
in comparison, much simpler beam generation hardware, 
namely one low-frequency phase modulator like a piezo- 
electric transducer, as opposed to two AOMs with RF 
driving electronics. On the other hand, the data pro- 
cessing for phase extraction is more complicated, which, 
however, becomes a smaller disadvantage with cheap pro- 
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FIG. 8. Angular resolution obtained by applying a DWS al- 
gorithm to the phases extracted from individual cells of a 
quadrant photodetector. 



cessing power. The heterodyne method typically requires 
additional stabilization loops [27l |28] to reach noise lev- 
els at pm/-\/Hz, e.g. for the laser power and certain 
common-mode pathlengths (see Ref. [S]). The exper- 
iments described above in Section IIVI show that these 
stabilizations are not required for the deep phase modu- 
lation technique. 



ments. This technique is based on the deep phase mod- 
ulation (over several radians) of one interferometer arm 
and can be considered as an extension of the well-known 
"Ji . . . J4" method [T51 - fT5] . The harmonic amplitudes are 
used to numerically solve an overdimensioned system of 
equations to extract the interferometer phase and other 
useful interferometer variables. This technique has been 
applied to experiments conducted on a very stable inter- 
ferometer (the engineering model of the LISA Pathfinder 
optical bench), achieving an optical pathlength readout 
sensitivity of the order of 20 pm/ VHz (which translates 
to lOpm/vHz for free-floating test mass displacement), 
and alignment measurements with an angular resolu- 
tion better than 10 nrad/vHz in the millihcrtz frequency 
band. This performance is comparable to the best hetero- 
dyne interferometers, and, e.g., only a factor of 2 above 
the LISA Pathfinder pathlength measurement require- 
ments. Two main noise sources were identified, namely 
laser frequency fluctuations and the frequency response of 
the analog portion of the data acquisition system, which 
both were completely mitigated by appropriate data pro- 
cessing methods, hence improving the performance of 
this technique by over a factor 35. Unlike other inter- 
ferometry techniques, no additional control loops, for in- 
stance, to actively stabilize the optical pathlength differ- 
ence or laser power fluctuations, have been implemented 
or are required to reach the current sensitivity. Nonethe- 
less, this could easily be done, in order to further improve 
the performance of this method. 



VIII. CONCLUSIONS 

We have presented an interferometry technique for 
high sensitivity length and angular optical measure- 
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